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STUB  PENETRATION  IN  CHEMICAL  MUNITION  SHIPMENT 

by 

Peter  S.  Westine 


Introduction 

The  M55  warhead  carries  a  chemical  agent  which  must  be  protected  from 
accidental  puncture  in  a  shipment  or  handling  mishap.  In  this  project, 
computer  runs  were  made  to  calculate  the  threshold  impact  velocity  if  a  pallet 
of  munitions  or  a  protective  package  containing  several  pallets  were  to  be 
dropped  onto  rigid  stubs  of  different  radii.  Figure  1  shows  schematically  the 
M55  pallet  with  15  rounds  each  inside  a  shipping/launching  tube.  The  shipping 
tubes  are  made  of  fiberglass;  therefore,  if  contacted  from  above  or  below, 
both  the  fiberglass  container  and  the  aluminum  walls  of  the  round  must  be 
penetrated  before  any  chemical  agent  leakage  can  occur.  Figure  2  presents  the 
cross  sections  which  must  be  penetrated  for  leakage  to  occur  under  these 
conditions,  which  shall  be  called  Case  1.  The  total  weight  of  one  pallet  with 
its  fifteen  rockets  is  1350  lb. 

Should  the  pallet  in  Figure  1  rotate  so  that  the  impact  is  on  the  side  of 
a  pallet,  extra  protection  is  provided  by  wooden  two-by-fours  in  the  pallet 
side  walls.  Thus,  the  Case  2  penetration  analysis  is  for  the  threshold  impact 
velocity  to  just  penetrate  a  composite  of  wood,  fiberglass,  and  aluminum  when 
the  pallet  is  dropped  onto  rigid  stubs  of  different  radii.  The  pallets  being 
dropped  in  Cases  1  and  2  are  identical.  Only  the  pallet  orientation  and 
number  of  layers  to  be  penetrated  change. 

The  third  accident  scenario.  Case  3,  is  puncture  of  a  shipping  container 
which  might  be  loaded  on  a  flat  bed  railroad  car  or  highway  vehicle.  These 
protective  packages  are  complex  containers  within  a  container.  The  one  shown 
schematically  in  Figure  3  and  analyzed  in  this  study  is  called  the  CAMPACT. 
Many  pallets  of  munitions  would  be  stored  in  a  CAMPACT  shipping  container; 
however,  in  Case  3  we  only  determine  the  threshold  velocity  for  just 
penetrating  a  CAMPACT  dropped  onto  different  radii  stubs.  Thus,  Case  3  does 
not  include  penetration  of  a  loaded  CAMPACT  container  plus  pallets  of 
munitions;  however,  as  a  first  approximation,  the  threshold  velocity  for  a 
pallet  inside  the  CAMPACT  is  the  square  root  of  the  sum  of  the  velocities 
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Figure  2.  Penetration  Conditions  in  Case  1 


Figure  3.  CAMPACT  Construction  Features 


squared  for  Case  1  (or  2)  and  Case  3.  A  CAMPACT  fully  loaded  with  eight 
pallets  weighs  approximately  46,100  lb. 

A  cross-section  through  the  walls  of  a  CAMPACT  is  seen  in  the  top-half  of 
Figure  4.  Many  of  the  layers  may  as  well  not  exist  for  puncture  analysis. 

The  steel  skins  surrounding  the  polyurethane  foam  are  basically  barriers  for 
containing  the  foam  when  it  is  foamed  in  place  after  the  inner  and  outer 
CAMPACT  containers  are  assembled.  The  outer  steel  sheet  is  actually 
corrugated,  but  it  is  only  12  mils  thick  and  is  used  for  keeping  out  rain. 
Kaowool  insulation,  a  very  weak  structural  material,  is  a  light,  thermal 
insulation  layer  installed  to  protect  the  container  contents  from  fire.  For 
all  practical  purposes,  the  major  anti-penetration  armor  is  the  Kevlar  as  our 
computations  establish.  The  polyurethane  foam  was  included  in  the  analysis 
because  il.  contributed  a  small  amount  of  resistance  due  to  its  large 
thickness.  The  actual  Case  3  cross-section  studied  computationally  is  that 
seen  in  the  bottom  half  of  Figure  4. 

In  this  report,  we  will  first  present  the  results  obtained  from  our 
calculations.  A  SwRI  penetation  code  under  development  for  the  U.S.  Army  Tank 
Automotive  Command  cal’ed  "Westine  Eight-Stage  Bodner  Type  Penetration  Model" 
was  used  in  all  calculations.  Details  of  how  this  code  functions  physically 
are  described  later  in  this  report. 

Results 

Figure  5  is  a  plot  of  the  critical  threshold  penetration  velocity  Vcr 
versus  stub  radius  R  for  Cases  1  and  2,  for  penetration  of  the  M55  chemical 
munition  in  a  pallet  as  in  Figure  1.  Failure  occurs  when  the  stub  punches 
through  the  aluminum  wall  so  that  agent  could  leak  out. 

Two  curves  are  seen  in  Figure  5.  The  lower  curve  is  for  Case  1 
penetration  as  sketched  in  Figure  2  without  any  wooden  crate  between  the 
shipping  tubes  and  the  stub.  The  upper  curve  is  for  the  Case  2  penetration, 
penetration  first  of  a  soft  pine  two-by-four,  followed  by  penetration  of  the 
fiberglass  shipping  tube  plus  aluminum  wall  of  the  munition.  Properties  for 
the  two-by-four  were  assumed  to  be  a  weight  density  of  0.0126  lb/in'*  (21.8 
lb/ft3),  a  flow  stress  of  1800  psi,  and  an  ultimate  strain  of  0.04  in/in. 

Woods  vary  greatly  in  material  properties  even  among  the  same  species.  Where 
the  wood  was  taken  from  in  a  tree,  whether  grains  are  close  together  or  far 
apart,  etc.,  all  influence  mechanical  properties.  These  particular  values 
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Figure  4.  Penetration  Conditions  for  Case  3 
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were  selected  as  being  representative  of  a  weak  pine.  A  weak  pine,  rather 
than  a  wood  with  stronger  properties,  was  used  so  the  calculation  would  be 
conservative.  In  all  of  these  computations,  the  stub  is  treated  as  long  and 
rigid. 

As  has  already  been  described,  the  third  case  is  penetation  of  a  CAMPACT 
shipping  container.  Figure  6  is  a  plot  of  threshold  penetration  velocity 
versus  stub  radius  for  threshold  of  penetration  into  the  46,100  lb  fully 
loaded  CAMPACT  container.  In  some  studies,  the  resistance  of  the  CAMPACT  to 
puncture  has  been  studied  by  remaking  the  container  into  an  equivalent 
thickness  of  mild  steel.  This  approach  is  not  one  which  we  favor  because  the 
equivalent  thickness  changes  with  stub  radius  as  well  as  mode  of  response. 
Nevertheless,  we  did  make  calculations  for  equivalent  threshold  mild  steel 
thickness. 

After  specific  values  for  CAMPACT  threshold  penetration  velocity  onto 
different  radii  stubs  had  been  obtained,  each  of  the  velocity  and  radius 
combinations  was  used  In  the  same  computer  program  to  determine  the  thickness 
of  rolled  homogeneous  mild  (1020)  steel  that  would  just  defeat  each  impact 
combination.  The  steel  was  assumed  to  be  1020  cold-rolled  with  a  flow  stress 
of  56,000  psi,  a  weight  density  of  0.283  lb/in^,  and  an  ultimate  strain  of 
0.15  in/in.  Also  plotted  in  Figure  6  Is  the  equivalent  thickness  of  a  1020 
mild  steel  armor.  These  thicknesses  change  with  stub  radii.  With  different 
physical  processes  taking  on  other  relative  importance  as  target  thickness  and 
material  properties  change,  use  of  the  equivalent  armor  concept  is  question¬ 
able  as  the  same  physical  processes  are  not  always  being  studied.  We  recom¬ 
mend  using  the  lower  curve  in  Figure  6  directly  wheneve**  system  assessments 
are  conducted. 

Computer  Code  Validation 

To  demonstrate  that  this  code  can  properly  handle  stub  penetration  into 
objects  at  low  velocities  of  impact,  two  types  of  comparisons  were  made.  In 
the  first,  we  went  to  the  literature  and  gathered  a  large  quantity  of  test 
data  obtained  by  Spaller*  in  1966.  In  one  series,  of  experiments  Spaller 
drop-tested  blocks  of  lead  onto  stubs  of  different  diameters.  Different 


*A.  E.  Spaller,  Structural  Analysis  of  Shipping  Casks,  Vol.  2  Resistance  to 
Puncture,  September  1966,  ORNL  TM  1312.  ~ 
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Figure  6.  Threshold  of  Penetrating  CAMPACT 


weights  of  lead  blocks  were  drop-tested  (either  26,  52,  78,  or  104  pound 
weights)  on  various  diameter  stubs  (either  3/8,  1/2,  or  5/8  inch  diameter). 
Velocity  of  impacts  ranged  from  6.7  to  15.4  ft/sec  for  52  pound  weights.  The 
recorded  response  for  different  impact  velocities  was  depth  of  stub 
penetration  into  the  lead.  Figure  7  shows  Spaller's  test  results  plotted 
nondimensionally  as  scaled  depth  of  penetration  in  stub  diameters  Y/d  versus 
scaled  energy  of  impact  E/ou-|td3.  By  scaling  the  test  results  all  can  be 
presented  on  a  single  graph.  The  ultimate  for  the  lead  was  17,550  psi  and  the 
weight  density  was  0.411  1b/1n3.  Because  lead  is  ductile,  the  flow  stress  is 
lower  closer  to  the  yield  point.  The  solid  line  in  Figure  7  was  calculated 
using  the  computer  code  used  in  this  study  with  a  flow  stress  of  6,800  psi  for 
lead.  Our  purpose  In  making  this  run  for  this  simplified  case  was  to 
demonstrate  that  fairly  accurate  results  can  be  obtained  using  this  analysis 
procedure  for  homogeneous  material  penetrations. 

The  second  comparisons  to  be  made  are  for  penetration  Into  Kevlar. 

Kevlar  is  made  by  bonding  together  fine  fibers.  Such  a  material  is  not 
homogeneous.  Parallel  and  perpendicular  to  the  fibers,  Kevlar  has  very 
different  tensile  as  well  as  compressive  ultimate  strengths.  Kevlar  29  fibers 
in  tension  have  an  ultimate  strength  of  approximately  400,000  psi.  In 
compression  parallel  to  the  fibers,  it  will  have  an  ultimate  strength  of 
approximately  76,000  psi,  and  perpendicular  to  the  fibers,  an  ultimate 
strength  of  approximately  38,000  psi.  Some  of  these  discrepancies  are  caused 
by  different  penetration  processes  occurring.  When  leaded  across  fibers,  some 
failures  are  associated  with  fiber  separation  and  spreading  rather  than  a 
tensile  tearing  of  fibers.  The  computer  code  that  we  used  is  currently 
designed  to  predict  the  response  of  homogeneous  not  heterogeneous  materials. 

A  representative  flow  stress  of  50,000  psi  was  used  for  Kevlar,  but  its  use 
needed  validation.  Large  amounts  of  unclassified  penetration  data  into  Kevlar 
is  unavailable;  however,  three  data  points  were  founded. 

The  designers  of  CAMPACT*  describe  a  one-quarter  scale  replica  model  test 
designed  to  meet  the  intent  of  10CFR7I.73.  Failure  energy  in  this  one-quarter 
scale  model  was  6,900  ft-lb.  Because  this  model  weighed  720  lb.  this  energy 
implies  that  the  threshold  velocity  of  impact  in  this  model  was  24.8  ft/sec. 

♦Information  communication  with  Robert  M.  Burgoyne,  6.  A.  Technologies,  Inc., 
December  1984. 


Our  computer  code  run  predicted  a  threshold  impact  velocity  of  25.1  ft/sec 
when  the  flow  stress  equaled  50,000  psi.  This  comparison  indicated  that  a 
reasonable  value  was  being  used  for  flow  stress,  but  one  point  comparisons  can 
be  fortuitous. 

Other  Kevlar  penetration  data  were  furnished  by  the  U.S.  Army  Mechanics 
and  Materials  Research  Center,  Watertown,  Massachusetts.*  Figures  8  and  9 
present  unclassified  ballistic  limit  data  for  5.56  mm  and  7.62  mm  ball 
ammunition  fired  Into  Kevlar  panels.  The  5.56  mm  ball  projectile  associated 
with  results  in  Figure  8  weighs  0.00766  lb.  At  an  initial  impact  velocity  of 
2400  fps,  we  calculated  the  ballistic  limit  thickness  to  be  0.683  inches. 

This  value  is  essentially  the  same  as  that  shown  in  Figure  8.  The  7.62  mm 
ball  projectile  associated  with  results  presented  in  Figure  9  weighs  0.02167 
lb.  At  a  muzzle  velocity  of  2740  fps  as  shown  In  Figure  9,  we  calculated  a 
ballistic  limit  thickness  of  1.38  inches.  This  too  is  essentially  the  same  as 
was  observed  empirically.  The  predictive  consistency  in  these  three  Kevlar 
comparisons  indicate  that  a  flow  stress  of  50,000  psi  is  a  realistic  value  for 
Kevlar  29. 

Computer  Code  Discussion 

The  computer  code  used  in  this  analysis  is  a  new  one  under  preparation 
by  SwRI  for  the  U.S.  Army  Tank -Automotive  Command,  Warren,  Michigan.  The 
approach  being  used  was  first  developed  by  S.  R.  Bodner  and  his  graduate 
students  at  the  Technion  in  Haifa,  Israel  In  1982  for  a  rigid  penetrator 
impacting  armor.**  Instead  of  writing  a  conservation  of  energy  equation  or  a 
force  equilibrium  equation,  they  wrote  an  energy  rate  equation.  The  rate  at 
which  kinetic  energy  leaves  the  penetrator  and  attached  effective  target  mass 
must  equal  the  rate  at  which  energy  is  dissipated  by  strain  energy,  frictional 
effects,  shear  along  boundaries  and  convective  inertial  effects.  As  will  be 
apparent  in  discussions  of  derivations,  this  rate  equation  Is  equivalent  to 
Newton's  equation  of  motion,  except  some  parameters  representing  the  extent  of 
a  dynamic  plastic  flow  field  must  be  obtained  by  differentiating  some  of  the 


*Massianica,  Frank,  Ballistic  Technology  of  Lightweight  Armors, 
AMMRC-Tr-76-15,  1976lconf  1dent1aT]k 

**M.  Ravid  and  S.  R.  Bodner,  "Dynamic  Perforation  of  Viscoplastic  Plates  by 
Rigid  Projectiles,"  Int.  J.  Enqnq.  Sc.,  Vol.  21,  No.  6,  pp.  577-591,  1983. 
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Figure  8.  Protection  Provided  by  Bonded  Kevlar  Fabric 
Against  5.56mm  Ball  Ml 93  Ammunition 
Army  Materials  and  Mechanics  Research  Center 


Figure  9.  Protection  Provided  by  Bonded  Kevlar  Fabric 
Against  7. 62nm  Bali  M80  Ammunition 
Army  Materials  and  Mechanics  Research  Center 


a  dynamic  plastic  flow  field  must  be  obtained  by  differentiating  some  of  the 
terms  to  min1ml2e  the  external  forces  applied  to  the  penetrator. 

In  addition,  a  second  key  earlier  contribution,  made  by  Bodner  and  other 
earlier  graduate  students*,  was  that  projectile  penetration  goes  through  a 
series  of  stages.  One-stage  engineering  penetration  models  are  seldom 
adequate.  In  our  analysis,  for  a  plugging  mode  of  failure,  we  use  a  six-stage 
penetration  process.  At  first  the  nose  must  become  embedded  In  the  target. 

In  thick  plates,  nose  embedment  and  a  subsequent  second  stage  result  in 
material  being  ejected  from  the  front  face  of  the  target.  The  top  two  figures 
in  Figure  10  illustrate  this  behavior.  Eventually,  as  the  depth  of  penetra¬ 
tion  grows  larger  (stage  3),  material  is  no  longer  ejected  from  any  target 
face,  and  the  penetration  process  becomes  a  cavity  expansion  process.  As  the 
penetration  process  increases,  the  plastic  flow  field  ahead  of  the  penetrator 
intersects  the  rear  face  of  the  armor,  and  a  bulge  begins  to  form.  The  stage 

4  bulge  grows  larger  until  the  ultimate  strain  is  reached,  and  a  plug  in  stage 

5  of  Figure  10  is  ejected.  Finally,  in  stage  6,  the  plug  is  ejected  and  the 
penetrator  begins  to  emerge. 

Whenever  the  armor  is  thin,  the  plastic  flow  field  ahead  of  the 
penetrator  will  Intersect  the  rear  surface  of  the  armor  before  the  nose  is 
fully  embedded.  Under  this  condition  a  stage  4  type  bulge  will  occur  which 
will  lead  to  petal ing  of  the  armor  when  the  ultimate  strain  is  exceeded  in 
stage  4.  As  soon  as  the  petal  begins  to  form,  the  armor  is  breeched.  This 
second  sequence  of  events  can  occur,  and  is  referred  to  as  a  petal ing  mode  in 
the  sequence  of  events  as  opposed  to  a  plugging  mode. 

To  understand  how  this  analysis  is  derived,  consider  a  rigid  penetrator 
in  a  stage  2  flow  field  as  shown  in  Figure  11.  In  Figure  11  a  plastic  flow 
field  extends  out  a  distance  NR  from  the  focal  point  for  the  nose  of  a 
penetrator  of  radius  R.  This  assumption  basically  states  that  a  plastic 
disturbance  extends  out  a  uniform  distance  of  (N-l)R  from  the  surface  of  a 
penetrator  with  a  hemispherical  nose.  Two  zones  are  seen  in  this  flow  field, 
zones  1  and  2.  Relative  to  the  target  itself  the  velocity  Vj _ j j  comes  from 
conservation  of  the  mass  flow  rate: 


*V.  Awerbuch  and  S.  R,  Bodner,  "Analysis  of  the  Mechanics  of  Perforation  of 
Projectiles  in  Metallic  Plates,"  Int.  J.  Solids  Structures,  Vol.  10,  pp. 
671-684,  1974. 


Figure  If).  Stages  in  long  Rod  Penetrator  Perforation  of  Thick  Target 
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Relative  to  the  penetrator,  V j- _ j j  is: 
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Now  a  self-consistent  flow  field  velocity  distribution  can  be  determined. 
Relative  to  the  penetrator,  the  radial  Vr  and  angular  V  velocities  are: 
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The  velocities  can  be  inspected  relative  to  the  flow  field  boundaries  to 
see  that  these  are  correct.  At  e  =  0°,  VQ  =  0  and  at  e  =  90°,  VQ  *  V j _ j j . 

At  e  =  0"  and  r  =  NR,  vr  =  -V;  at  9  =0°  and  r  =  R,  Vr  =  0;  and  at  e  =  90°, 

Vr  *  0.  In  zone  II,  the  motion  is  rigid  body  movement  with  the  top  and  bottom 
moving  at  the  same  velocity.  All  velocities  are  relative  to  a  coordinate 
system  attached  to  the  penetrator. 

Now  that  the  velocity  field  has  been  obtained,  everything  follows 
algebratically.  The  first  energy  rate  to  be  calculated  is  the  strain  energy 
rate  wv.  To  do  this  the  strain  rate  fields  must  be  obtained  from  the  velocity 
field  by  differentiation.  The  strain  rates  in  zone  I  are: 
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In  zone  II,  rigid  body  motion  results  so  that 
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The  differential  volume  in  zone  I  i s : 
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Strain  energy  rate  wy  in  the  plastic  flow  field  is  given  by: 
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which  is: 
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After  substituting  Equations  (5)  and  (7)  into  equation  (9)  and  integrating  one 
obtains: 
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There  is  no  w,.  in  zone  II  because  the  e's  are  zero. 
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Around  the  nose  and  over  the  embedded  length  of  the  penetrator  there  is  a 
frictional  shear  energy  rate  w^.  This  wf  is  given  by: 
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where  Avi  is  the  velocity  of  one  surface  relative  to  the  other  and  dA  is  the 
differential  area  of  the  interface.  If  t  Is  assumed  to  be  no,  where  u  is  a 
coefficient  of  friction,  the  friction  energy  rate  is: 

2  M2  .  M2  v 

wf  S  J  U0  — 5 -  V  sin  6  2 n  R  sin  6  de  +  po  ■  2tt  R(Z-R)  (12) 

o  (N  -l)  (N  -1) 


or 


2  iro  V  R2[  y  (Z/R  -  0.2146)  — 5 -  ] 

(N-1) 


(13) 


In  addition,  a  shear  energy  rate  w  occurs  along  boundaries  of  zones 

s 

because  they  move  relative  to  one  another.  The  quantity  w  is  given  by  the 

s 

same  general  equation  as  wf,  equation  (11).  In  stage  2,  ws  is  given  by: 
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The  shear  stress  the  Huber-Henke-Mises  criteria,  relates  shear  stress  to 
normal  yield  stress.  Completing  the  integration  gives: 
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Before  inertial  energy  rates  can  be  determined,  the  velocity  equations, 
equation  (4),  must  be  differentiated  with  respect  to  time  to  obtain 
accelerations.  In  zone  I,  the  accelerations  are 
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Notice  that^two  types  of  acceleration  terms  exist:  those  with  V  in  them 

and  those  with  —  in  them.  This  means  that  two  different  types  of  inertial 
r 

energy  rate  terms  will  exist.  The  convective-inertial  terms  w  are  obtained 

V  * 

by  using  only  those  acceleration  terms  with  —  in  equation  (16).  The  dynamic 
acceleration  terms  are  obtained  by  using  only  the  v  acceleration  terms  in 
equation  (16).  No  matter  which  type  of  inertial  energy  rate  is  being 
computed,  the  inertial  energy  rate  is  given  by: 
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Substituting  the  velocity  and  acceleration  fields  into  equation  (17)  gives  for 
the  convective  energy  rate  w^: 
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The  dynamic  acceleration  energy  rate  wd  is: 
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The  final  energy  rate  term  is  the  energy  rate  in  the  rigid  penetrator. 


The  penetrator ‘s  inertial  energy  rate  term  wp  is: 
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Now  all  the  energy  rate  terms  have  been  derived  for  stage  2  in  the  pene¬ 
tration  process.  Notice  that  of  all  the  energy  rate  terms,  only  equations 
(21)  and  (22)  have  V  V  in  them.  Equation  (21)  is  the  local  dynamic  accelera¬ 
tion  which  adds  an  effective  mass  of  target  material  to  the  mass  of  the  pene¬ 
trator.  If  energy  rates  are  equated,  the  energy  rate  loss  from  the  sum  of 
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equations  (21)  and  (22)  must  go  Into  wg  +  wf  +  w  .  Mathematically  this 
process  Is: 


[H  .  2,  p  R3{!iV-6,N2  ,<■  Sg-^aNV-O)]  . 

9(n -ir 


.!».  I4l«ir  -'2pj,-n-N  o  R2  -  2*  y(|  -  0.2146)  — jp—  a  R2 
(m2-i)  r  (n2-i) 


-[0.4534  —3 -  +  0.5773  — 5 -  (|  -  1)]  2-ir  0  R2 

(N2-1)  (N2-1)  R 


r  it  N2f2N2-1 ) 

^  (n2-i)2 


R2  V2 


(23) 


By  substituting  v  for  ^  .  equation  (23)  could  be  solved  for  dV  using 
incremental  dz  steps  if  the  extent  of  the  plastic  flow  field  given  by  N  could 
be  obtained.  Equation  (23),  obtained  by  balancing  energy  rates  is  Newton's 
equation  of  motion.  The  left-hand  side  is  the  total  acting  mass  (penetrator 
and  effective  target  mass)  times  the  acceleration.  All  of  the 
wv  +  wf  +  wg  +  terms  on  the  right-hand  side  give  the  external  forces 
acting  on  the  penetrator.  This  point  is  a  very  important  statement. 

To  obtain  N  so  that  the  equation  of  motion  can  be  solved,  we  differen¬ 
tiate  the  right-hand  side  of  equation  (23)  only  with  respect  to  N  and  set  the 
result  equal  to  zero.  This  process  says  that  the  extent  of  the  plastic  flow 
field  will  be  such  that  the  transient  external  forces  on  the  projectile  are 
minimized.  In  other  words,  the  projectile  follows  the  path  of  least  resis¬ 
tance.  Differentiation  thus  gives  a  second  relationship. 
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For  specific  values  of  Z/R  and  ^!L-  ,  equation  (24)  can  be  solved  for 
N.  After  N  Is  obtained  dV  Is  obtained  from  equation  (23).  Solving  equations 
(23)  and  (24)  are  Iterative  procedures.  Equation  (24)  Is  for  the  extent  of 
the  plastic  flow  field;  equation  (23)  Is  the  equation  of  motion. 

Similar  mathematical  procedures  were  followed  for  all  of  the  stages.  If 
the  equations  for  all  stages  were  presented,  more  than  200  equations  would  be 
presented.  By  presenting  the  mathematics  In  stage  2,  the  reader  sees  the 
processes  which  are  considered.  Those  wishing  complete  details  and 
derivations  for  all  stages  should  obtain  SwRI's  report  for  TACOM  when  it 
appears.*' 

All  threshold  penetration  velocities  were  obtained  using  the  same 
procedure.  If  a  target  was  multi-layered,  as  with  the  aluminum  munition  in 
the  fiberglass  shipping  container,  the  residual  velocity  of  the  stub  exiting 
the  fiberglass  became  the  input  velocity  into  the  aluminum.  Various  impact 
velocities  were  tried  in  a  trial  and  error  procedure  until  the  velocity  for 
just  penetrating  the  aluminum  was  within  0.1  ft/sec  of  a  velocity  which  does 
not  penetrate.  Finally  the  threshold  velocity  was  obtained  by  averaging  these 
two  velocities. 

Whenever  an  equivalent  thickness  of  steel  was  obtained,  various  mild 
steel  thicknesses  of  armor  were  tried  for  that  radius  of  stub  and  impact 
velocity  until  a  thickness  was  obtained  for  penetrating  and  another  for  not 
penetrating  which  were  within  0.01  inches  of  each  other.  Hence,  all  plots  of 
stub  radius  versus  impact  velocity  or  equivalent  thickness  of  mild  steel 
versus  stub  radius  come  from  a  sequence  of  trial  and  error  runs. 

The  advantage  of  analysis  procedure  relative  to  hydrocode  is  that  only 
two  equations  must  be  solved  at  each  time  step  rather  than  hundreds,  or  even 
thousands,  of  equations  as  in  any  numerical  hydrocode.  The  number  of 
equations  ic  decreased  by  using  this  engineering  analysis  technique. 
Nevertheless,  essentially  all  physical  processes  except  strain  rate  effects 
nad  penetrator  erosion  are  included.  That  the  analysis  can  work  is  already 
demonstrated  in  Figure  7  for  drop  tests  of  lead  blocks  onto  stubs. 

A  listing  for  the  entire  computer  program  will  appear  in  SwRI's  report 
for  TACOM  when  it  is  approved  for  release  by  the  sponsor. 


*P.  Westine,  Development  of  an  Engineering  Computer  Code  for  Eroding 
Projectile  Penetration  and  Perforation  of  Armor,  Contract  DAA07-84-C-R120. 
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